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SUMMARY 


Fundamental equations of conservation of mass, momentum, and 
energy are applied to the problem of heat transmission by forced 
convection from a line combustion source. The influence of atmos- 
pheric stability and surface roughness on growth of the heated zone, 
decay of temperature downstream from the heat source, and aerodynamic 
drag is presented. Theory is general, being based on the eddy diffu- 
sivity for momentum given by Deacon (2)- which is 



Eddy diffusivity for heat is taken as equal to that for momentum. 

Theory presents groupings of the variables that eire common to 
heat transfer and fluid flow problems concerning the surface air layer. 
Results reveal that atmospheric stability plays an important role in 
the growth of the heated region, decay of temperatirre, and aerodynamic 
drag for small values of friction velocity and macroviscosity parameter. 
For large values of these parameters (large surface roughness), atmos- 
pheric stability is not important. 


T7 For underlined niambers in parentheses, refer to literature 
cited, page 40. 
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INTRODUCTION 


Combustion of a fuel generates heat which must be transferred. 
Transmission of thermal energy is by the conventional modes of heat 
transfer — conduction, free and forced convection, and radiation. A 
line fire burning with a moderate wind speed in light fuels, such as 
grass, the heat transfer is predominately by forced convection. On 
such a fire as well as on some back fires, free convection plays a 
minor role since the momentum transfer by the buoyant forces is small 
compared to the momentum transfer by the inertia forces existing in 
the velocity field. Conduction to the groimd is negligible. Flame 
radiation which depends upon flame characteristics and geometry is of 
greater importance. Experiments were conducted by the Forest Products 
Research Laboratory in England (9) on the combustion of crib fires with 
wood as a fuel. Heat balances made on the combustion process indicated 
that heat transferred by convection was 2 to 4 times greater than by 
radiation. If the flame geometry is small, in 5 )lying moderate rates of 
energy release, the radiation exchange may be neglected. On the basis 
of the foregoing postulates, the convection zone may be studied by 
using heat transfer and fluid flow theories. 

In a study of the convection zone generated by a combustion 
process, the combustion mechanism may be neglected. Heat generated by 
combustion of the fuel is assumed to be known and from the theoretical 
viewpoint represents the heat source strength. Theory is based upon 
conservation of mass, momentum, and energy which is amenable to mathe- 
matical analysis. Despite its complexity, the combustion process is 
itself of theoretical Interest. To understand a fire burning in the 
atmosphere, the steady-state combustion mechanism must be studied. 
Ultimately, the transient combustion process requires investigation. 
Since the power of any combustion process, whether it be for construc- 
tive or destructive purposes, is determined at the combustion zone, 
this region is of prime importance. It is not the purpose of this 
report to discuss combustion theory but only to stress its importance 
and direct relation to heat transfer. The combustion process generates 
sensible heat which must undergo transmission. By assuming the com- 
bustion zone as a line source of heat, transfer of heat by forced 
convection may be analyzed. A generalized approach is presented inves- 
tigating forced convection heat transfer over rough surfaces as 
influenced by atmospheric stability. 

Prior to analysis of the specific problem, the fundamental equa- 
tions are presented including the buoyancy term in the momentum equation. 
For a given heat source strength, the minimum magnitude of wind velocity 
is derived so that the buoyancy term may be neglected. When this condi- 
tion is imposed a solution of the system of equations is possible. 

Thereafter the system of partial differential equations is 
transformed to a system of ordinary differential equations. Solutions 
with appropriate boundary conditions are obtained using the differential 
analyzer. Growth of the heated zone may be determined as a function of 
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friction velocity ratio, macroviscosity parameter, Reynolds Number, 
and atmospheric stability. Temperature profiles and decay of profiles 
downstream from the heat source as a function of velocity ratio, macro- 
viscosity pEirameter, Reynolds Number, and atmospheric stability are 
presented. The heat source is expressed in terms of the combustion 
process. Finally, the shear stress and drag coefficient are determined. 

Results of the investigation are applicable to line fires of 
moderate heat source strength wherein forced convection is the dominat- 
ing mode of transfer. Likewise, the results are of interest in the 
diffusion of mass from a line source which is of importance in chemical 
warfare, smoke screen operations, and disposal of industrial plant 
waste gases. 


ANALYSIS 

Basic Equations 

The continuity, momentum, and energy equations for the two- 
dimensional flow due to a line source for a fluid of constant Prandtl 
Number, specific heat, and viscosity are 


hn hw_ 


Su Su Sp 


Sw Sw 


+ ^ 


PgJCpU ^ + pgJc^w 


aT 
p“'5S 



(1) 

a^u a^u'\ 
ax^ az"/ 

(2) 

a^w a^w \ 

. ax^ "" az^ / 

(3) 

Ux^ az^/ 

(4) 




In most meteorological problems the pressure change in the 
direction of the mean wind is small, especially within the surface 
boundary layer region. This term is small compared to the other 
terms and may be neglected. 


Applying equation (3) in the region exterior to the zone of 
influence states 




The foregoing equations become 


Su 



(5) 
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From the perfect gas law at constant pressure 


= 1 . 
p " T. 


there results 


P„ - P = P 
00 


T - T \ 
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00 

The foregoing equations become 
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' 00 oo 
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(9) 

( 10 ) 

( 11 ) 

( 12 ) 


In addition, the relation stating that the heat output of the 
source is accounted for at any cross-section in the region of influence 
yields 


I -/ 


(13) 


The theory may now be expressed in nondimens ional form. If 

is a typical velocity and x a typical length, then — — - becomes 

00 

a nondimens ional parameter. Writing 
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X = — 

X 
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The boundary conditions are 



z = 0 

u = w = 0 

i=° 

Sg , 
Sz 
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00 
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00 
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Upon substitution the equations become 
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On the assumption that the derivatives which occur remain finite 

X 

when — » 1, the limiting form of these equations is 


o 

II 

+ 

(14) 

„ ,, au 

(15) 

SP . Q 

^ Lp c T 

00 p 00 00 

(16) 

,, aa ^ aG i a^a 

(IT) 

p CO 

/ UGdZ = 1 

^ n 

( 18 ) 
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From equation (16), restriction to a boundary- layer type problem 
which neglects the buoyancy term is given by 

^ « p c T u^ 

L 00 p 00 00 

The function G is of order unity, having a value of zero at the 
exterior region and 0.66 at the surface from equation (34) • Property 
values of air to be used throughout this paper are 


1.80 X 10""^ 

sq.ft/sec . 

2.34 X 10"^ 

lb. sec.?' ft.'^ 

0.24 BTU/lb. 

°R 

530°R 



Using the foregoing values the relation is 

^ « 0.30 u^ 

L 00 

When the term is arbitrarily reduced_to 10 percent of the pre- 
ceding value, that is, letting Q/L = 0.03 u , the relation is shown in 
figure 1. For a known heat output, the velocity must be equal to or 
greater than the value given by figure 1. Buoyancy may then be neglected. 
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Figure 1. — Minimum velocity for given heat output to 
neglect buoyancy 


The relation has been calculated for a line combustion of a forest fuel 
of rate R tons per acre per minute and width of flame front L*. 

This is shown in figure 2. 



Figure 2. --Minimum wind velocity for given fuel 
consumption to neglect buoyancy 
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The meaning of the relationship shown in Figure 1 (or the relation- 
ship in Figure 2) can he illustrated by using it to compute specific rates 
of spread corresponding to a given wind velocity for which bu^ancy can 
be neglected. For example, in a cover type with a light fuel— of 0.5 
tons per acre in which a fire is spreading at the rate of 30 feet per 
minute, the wind speed would have to be 14 feet per second or greater in 
order that buoyancy could be neglected. For this same light fuel, 0.5 
tons per acre, but with a rate of spread of 3 feet per minute, the wind 
speed would have to be 6.5 feet per second or greater for buoyancy to be 
neglected. However, if the wind speed were 6.5 feet per second for a 
considerably heavier fuel, 15 tons per acre, then the rate of spread of 
a fire in this fuel would have to be 0.1 feet per minute or less for 
buoyancy to be neglected. 


The condition for neglecting buoyant forces has been established 
in the preceding section. To apply equations (9) to (12) inclusive to 
a transport process in an atmosphere it is necessary to introduce vari- 
able diffusivities of momentum and heat and Kjj(z). The equations 

(9), (10), and (12) in dimensional form then become 


^u , Sw 


= 0 


( 


u 


Bu 


+ w 


Su 


d 


K (z) 
m' ' 


Su 


a(T - T ) 


S(T - T ) 


u 


dz 


+ w 


dz 


a 


yz) 


a(T - T )■ 


(19) 

( 20 ) 


( 21 ) 


with boundary conditions 

S(T - T ) 

z = 0 u = w = 0 ^ — — = 0 


z — ><» u = u 
00 


T - T = 0 
00 


Deacon (2) studied the vertical diffusion in the lowest layers of 
the atmosphere. Field experiments were conducted over two types of natural 
land surfaces. From the meas\irements, the vertical component of the eddy 
viscosity was determined for the lowest 30 meters of the atmosphere. The 
result was 


K (z) = k V 
m 00 



( 22 ) 


27 Computation based on a heat value of fuel of 8,000 BTU per 

pound. 
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Rider 0) conducted an experimental study of the factors which 
influence the vertical turbulent transport of momentum, water vapor, 
and heat in the first two meters of air above a short grass surface. 
Measurements were made of the vertical profiles of velocity, humidity, 
and temperature, evaporation rate, and aerodynamic drag of the surface. 
An attempt was made to measure the individual terms involved in the 
heat balance. 


Results indicated that the wind velocity and temperat\rre profiles 
always have the same distribution, and this distribution is representa- 
tive of the humidity profile in the majority of cases. The surface 
friction coefficient was found to be independent of wind speed, increas- 
ing in magnitude during unstable conditions. Values of the eddy dlffu- 
sivity for momentum, water vapor, and heat were obtained by measuring 
the quantities in the expressions defining the diffusivities . Results 
of this investigation showed that the diffusivities for momentum and 
vapor are identical over the range of stability of the experimental 
investigation. The value of the diffusivity for heat was approximately 
the same as for momentum and vapor, except that on occasion it is much 
larger than the other two. 


Despite limited experimental evidence, expressions for the eddy 
diffusivities for momentum and heat were taken as equal. Thus 


K„(2) = k 



( 23 ) 


The set of partial differential equations can be greatly simpli- 
fied by transformation to a new coordinate system. Introduction of the 
stream function t automatically satisfies the continuity equation. 

The velocity components are given by 


u 


Silf 



( 24 ) 


The formulas for making the transformations from the independent 
variables x,z to an arbitrary and completely new set x,t) are 


u 

n = — k 
‘ V 

00 



3 



2-3 

/ ^"*^0 



V ^00 . 

00 


3 


il§ 


/ U..Z 

2-3 
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2-P 


u X 
; 00 


1 

2-p 


( 25 ) 


ilf = k 


V , — , — 


U X 


2-D 


V ^00 / \ 


f(r\) 
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After transformation 


u = u 

on 


Sf 

00 OT) 


ll§ 


. -§_ .;h&] _ _L 

1 2 -B V2-6/ 2 -B 

/ u„ \ / u„ z \ - u X \ 

.(n) 


T - T = k 
00 


2-3 Q 


Lp gc V I u , 

00^ p 00 \ 00 / 


and equations (20) and (21) become 


/.P 


d i 




^ 1 , „ d"f 


dn^ I2-3; ^ ^^2 


= 0 


(27a) 


01/ W-i-' ilM = 0 

dn v2-3y dn 

with boundary conditions 

n=0 f=f'=o 

n — > oo f ' = 1 


(27b) 


g' = 0 
g = 0 


Conservation of energy at any cross-section downstream from the 
heat source yields 


n 00 

J Poog^p"^ (T - T^) dz = Q/L 


( 28 ) 


which is transformed to 


df ^ T 

15 ; 2 1 


( 29 ) 


Laminar Solution 


The laminar flow case is a special case of the more generalized 
system and is presented here in the interest of completeness. If 3=0 
/ u^z ■) 

and k [ / equals unity, then the eddy viscosity relation is that for 

V / 
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laminar flow (ct = l) 


y z) = y z) 


The equations 

3 2 

d^f 1 „ d^f 

dT) drj 


£s ^ i^fal 
2 d, 


= V 

00 

( 27 a) and ( 27 b) become 
= 0 

= 0 


(30a) 

(30b) 


Equation (30a) has been solved for the boundary conditions of 
this problem and is presented in Table III of reference (5). Solution 
of equation ( 30 b) is 


g(Tl) 



where K is given by equation (29) as 


( 31 ) 



0.66 


(32) 


By selecting the edge of the bovmdary layer as the position 
where the velocity attains 99 percent of the value of u / there 
results Tig =5.0 which gives “ 


6 u 
00 

V 

00 


/u x\ 1/2 

\ V ! 

00 / 


(33) 


The temperature at the surface downstream from the heat source is 


w 


T = 0.66 


Q 


■ u 


1 

2 


Lp gc V I 

'^00” p 00 \ 


(34) 


The laminar flow result states that the temperature is linearly 
dependent upon the heat source strength and decays inversely as the 
square root of the velocity and distance from the source. 
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Turbulent Solution 


Transfer processes in the atmosphere are not governed by laminar 
flow theory. In the foundation of the theory, the diffusion mechanism 
applicable to the atmosphere was used. 

Solution of the general system, equations (27a) and (27b), was 
undertaken by the University of California at Los Angeles, Department 
of Engineering Differential Analyzer.—' 


Integration of equation (27a) gives 




2 

dq^ 



f 


dX^ dq^ 


q=e 


( 35 ) 


Integration is started a small distance e from the origin due 
to the singularity at q = 0 . The value of the right side of equation 
(35) is defined as 


2 

.P jO 

dq 


= K, 


q=6 


( 36 ) 


In order to eliminate the two point boundary condition thereby 
facilitating machine solution the transformation 

F = f/A and N = q/B 

is introduced into equation ( 35 ) and choosing 
A = ( 2 -p)B^“^ and K, = AB^"^ 


there results on integration 



with bovindary conditions 



N = e/B 


I 

F and F finite 


dN 


N — > 00 


B/A 


( 37 ) 


37 Messrs. Malcolm R. Davis and Don Lebell are responsible for 
the mathematical steps which lead to setting up the equations for the 
differential analyzer including the trial and error method. 
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To 'begin the machine solution, the function f is expanded 
in a Taylor series about the point e . 


I » I 


f (ti) = f (e) + f (e) (n-e) + ^ + . . . . 

n+1^ 
nl 


+ (n-e) + . . . 

which for t) = 0 yields 


f'(0) = 0 = f'(e) - f"(e)e + " + 


f"'(e) 2 


+ (- 1 ) 


n n 


ni 


G 


or 


0 = 


f (g) 




,n n-1 


+ (-1) T-^ e + 


ni 


If the higher order derivatives are assumed finite and g is 
selected arbitrarily small, then the higher order terms may be neglected 
so that 


= ( 38 ) 

A value of g/B is arbitrarily selected with zero initial values 
of F and F . A run with the differential analyzer is completed 
which determ|i;ies A, B, and K,. From these results g may be calcu- 
lated and f (g) is determined by equation ( 36 ). Successively smaller 
values of e/^, are assigned and machine solutions obtained until}, the 
calculated f (g) when ins9rted in equation (38) yields an f (g) 
which is small compared to,f («) = 1. The final assumed value of g/B 
and calculated value of (f (g) B/A) with F as zero is run to see that 
the behavior of F is reproduced. 

The equation (27b) is transformed and integrated as 
fN F 

In G = - / dN ( 39 ) 

^ g/B 


where 
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The differential analyzer run was performed for three values of 
P , namely, 3/4^ and 5/4 which represent subadiahatic, neutral, and 
superadiabatic conditions . 

The condition of conservation of energy given by equation (29) 
is transformed to 

g(0) 

aJ § aw 

df* 

Results of the integration are the functions f and which 

are presented in figures 3 and 4* The function G is presented in 
figure 5. Also revealed in figure 6 is the integrand which appears in 
the conservation of energy relation. 


RESULTS 


Velocity Ratio and Macroviscosity Parameter 

Theoretical investigation of heat transfer from a line source in 
the atmosphere yields the effect of friction velocity, macroviscosity 
parameter, and atmospheric stability on the growth of the ten 5 )erature 
zone. For a given surface and magnitude of wind velocity the effect 
of atmospheric stability on the heat transfer may be determined. This 
effect is due to the variation of the vertical component of eddy diffu- 
sivity as the stability changes. Results are presented using values of 
friction velocity and macroviscosity given by Sheppard (6) . The friction 
velocity is defined by the relation 



Near the surface the shearing stress is considered constant and equal 
to the value at the surface. The friction velocity is then 



and depends on the nature of the surface and the magnitude of the mean 
velocity. The macroviscosity parameter is considered as a quantity 
which plays a role in fully rough flow analogous to that experienced 
by the kinematic viscosity in smooth flow. 
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^ GRADIENT OF FUNCTION (DIMENSIONLESS) f FUNCTION (DIMENSIONLESS) 



Figxire 3 •--Variation of function with independent 
variable as dependent on stability 



17 INDEPENDENT VARIABLE (DIMENSIONLESS) 

Figure 4 . — Variation of gradient with independent 
variable as dependent on stability 
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~ FUNCTION (DIMENSIONLESS) 6 FUNCTION (DIMENSIONLESS) 



77 INDEPENDENT VARIABLE (DIMENSIONLESS) 


Figure 5. 


— Variation of function with independent 
variable as dependent on stability 



Figure 6. — Variation of function with independent 
variable as dependent on stability 
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Growth of Heated Zone 


Growth of the heated zone is given by the generalized equation 

P 1-B 1 


u 6 
00 


2-p 


-7- = ^ k- 




2-p 


u 




izP 

2 -p 2 -B 

/U x^ ^ 

a 


( 40 ) 


The parameter p is dependent on the stability of the atmosphere and has 
been correlated with the Richardson Number by Deacon ( 2 ) . Figure 7 is 
the correlation taken from Deacon. Values of ti are determined from 
solution of the differential equation for the three values of p . 

These values are presented in table 1 . 


Table 1 . --Variation of Tig with p 


p 

^6 

0.75 

7.10 

1.00 

6.67 

1.25 

5.82 


Application of the theory requires specification of the Richardson Number. 
The stability parameter p is then obtained from figure 7. Figure 8 
yields the required magnitude of . For the natural surface considered, 
the value of the dimensionless friction velocity ratio and the 

macroviscosity parameter u^z ,v are available in the literature or must 
be determined experimentally. 'Growth of the region of influence 6 may 
be calculated as a function of distance from the source x for the speci- 
fied wind velocity u . 

00 

Sheppard (6) published values of friction velocity ratio and macro- 
viscosity parameter. These are presented in table 2 . 


Table 2 . --Representative values of uVu and u„z /v for 

00 * O «> 


natural surfaces. Neutral stability u^ corres- 

ponds to u = 500 cm/sec at 200 cm 

00 





00 

^*^o/v 
' 00 

Very smooth (mud flats, ice) 

0.0316 

0.096 

Lawn, grass up to 1 cm high 

0.050 

15.6 

Downland, thin grass up to 10 cm high 

0.071 

156.0 

Thick grass, up to 10 cm high 

0.089 

623.0 

Thin grass, up to 50 cm high 

0.107 

1647.0 

Thick grass, up to 50 cm high 

0.127 

3353.0 
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/9 stability parameter (DIMENSIONLESSl 


Ri RICHARDSON^ NUMBER (DIMENSIONLESS) 
-0.4 -0.2 -Ql 0 006 008 O.IO 



Figure 7. --Variation of p with stability 
(after Deacon) 



ST/«ILITY parameter (DIMENSIONLESS) 

Figure 8. --Variation of functions 

and g(0) with stability” 
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Theoretical results are presented in figure 9 which reveals 
growth of the heated zone. Natural surfaces of lawn grass up to 1 cm 
high and thick grass up to 50 cm high are compared. ‘A wind velocity of 
17 ft/sec was arbitrarily selected. For these conditions the effect 
of stability on growth of the heated zone is evident. The generalized 
correlation for the two surfaces is presented in figure 10. 

Temperature Decay 


Decay of temperature downstream from the heat source is given 
by the generalized equation 


T - 


T 

00 


g(0) 

1 




(41) 


The nondimens ional velocity ratio and macroviscosity parameter also appear 
in this equation. Strength of the heat source must be specified. The 
function G(ti) results from the solution of the differential equation 
and represents the generalized temperature profile in the heated layer. 

The surface temperature is maximum and is denoted by G(t)) = 1 . Values 
of g(0) for the three cases of p are presented in table 3 . 

Table 3. --Variation of g(0) with g 


3 

g(o) 

0.75 

0.70 

1.00 

0.94 

1.25 

1.58 


This fimctlon is also shown in flginre 8. Specification of the parameter 
gives the surface temperatxire as a function of distance downstream from 
the heated source for a specified wind velocity. Influence of the para- 
meters on decay of surface temperature is presented in figure 11. 

Natijral surfaces of lawn grass up to 1 cm high and thick grass up to 
50 cm high were selected. A wind velocity of I7 ft/sec and heat output 
of 60 BTU/sec ft were arbitrarily selected. This heat output corres- 
ponds to a forest fuel rate of 10 ton ft/acre min, based upon a 100^ 
combustion efficiency considering forced convection only. This fuel 
rate is 21.6 x 10^ BTU/hr ft and figure 1 indicates that a wind velocity 
of 13 ft/sec or greater is required permitting the neglect of the buoy- 
ancy effect. The influence of stability on decay of surface temperature 
is calculated. Generalized correlations for these two natural surfaces 
are presented in figure 12. 
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8 HEIGHT OF ACTION ZONE (FEET) 



0 100 200 300 400 


X DISTANCE DOWNSTREAM FROM SOURCE (FEET) 

Figure 9. — Grovth of action zone as function of velocity ratio, 
macroviscosity parameter, and stability parameter 


HEIGHT MODULUS (10 UNITS) 



Figiire 10. — Height of heated zone as function of velocity ratio, 
macroviscosity parameter, and stability parameter 
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surface temperature difference (®f) 



Figure 11. — Decay of surface tenqperature as function of velocity ratio, 
macroviscosity parameter, and stability parameter 



Va, 

Figure 12. — Decay of temperature as function of velocity ratio, 
macroviscosity parameter, and stability parameter 
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Decay of temperature downstream from the heat source may be calcu- 


lated in terms of the fuel evaporation rate M which represents the 
actual rate at which fuel is consumed. This rate does not represent 
radiative or convective energy rates directly since the efficiency of 
combustion and partition of energy are unknown. By a heat balance 
there results 


where L is the width i of the flame front, e is the combustion effi- 
ciency for forced convection, and H is the heat of combustion of the 
fuel. Equation (4I) becomes 



This expression relates the heat transfer to the combustion process. The 
fuel evaporation rate (weight rate per unit area) and the percentage of 
this weight rate that is converted into sensible heat for transfer by 
forced convection are required. 

Shear Stress and Drag Coefficient 


Defining the shear stress as 



( 43 ) 


yields the drag coefficient as 



T 



00 


o 


z — > o 


( 44 ) 


applying the transformation there results 



K 


( 45 ) 


- 29 - 


Using the definition of the friction velocity, the expression may be 
presented in general form as 


T 

C — ^ 

1 2 
^ u 

(C 00 00 




u 

. 00 - 


4P-4 

2-p 


* o 


2-2p 

2-p 



(46) 


Values of K, for the three selected cases of p are given in table 4 
and presented in figure 13. 


Table 4* — Variation of K, with p 


p 

K, 

0.75 

0.242 

1.00 

0.107 

1.25 

0.047 


For the case of neutral stability (p = l), the drag coefficient is 


T 

n _ ^ 

1 2 
oP u- 

d 00 00 


0.0036 


(47) 


This result states that the drag coefficient is independent of wind 
velocity, position, and type of natural surface. For cases of stability 
not equal to unity, the drag coefficient is dependent on wind velocity, 
position, and type of natural siirface. The shear stress at the surface 
is 


8/5 

T U ' 

O 00 


T = 0.0018 p u^ 
O 00 00 


T 

O 


8/3 

U • 


p = 0.75 
p = 1.00 
p = 1.25 


(48) 


Theoretical results of the shear stress are presented in figure 14 for 
natural surfaces of lawn grass up to 1 cm high and thick grass up to 50 cm 
high. The downstream position for the cases of stability unequal to unity 
was arbitrarily selected as 100 ft. Experimental values of the shear 
stress as determined by Vehrencamp (8) are indicated. 
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Figure 13 .--Variation of function K, with stability 



Figure 14. --Shear stress as a function of velocity 
(p = 1 independent of velocity ratio 
and macrovisGOsity parameter) 


DISCUSSION 


A mathematical theory has investigated forced convection heat 
transfer from a combustion process in the atmosphere. From this theory 
no discussion can be made of the influence of the atmospheric variables 
on the combustion rate, that is the buildup or decay of a fire. A 
discussion can be presented concerning the factors which influence the 
transfer of heat from a combustion zone considering only steady state 
forced convection. 


The initial result of the investigation is the condition which 
must be fulfilled so that the buoyancy effect may be neglected. Physi- 
cally, this is an energy relationship and is expressed as 


Lc T 
P ' 


« p u^ 
00 00 


This condition states that the energy rate per unit area of the combus- 
tion process must be much less than the kinetic energy rate per unit 
area of the existing velocity field. For a given combustion zone the 
wind velocity must be greater than a specified magnitude so that the 
heat generated will be transferred by forced convection. The case of 
free convection only has been treated in reference ( 4 ) . 


Theory was developed considering three cases of atmospheric 
stability. Strong subadiabatic, neutral, and strong superadiabatlc 
conditions were selected. Variation of the stability parameter taken 
from Deacon (2) is shown in figure 7- Values of p of 0,75 and 1.25 
represent extreme cases. 


The velocity ratio and macroviscosity parameter are determined 
by the type of natural surface. Further discussion of these quantities 
is presented in references (6,7>8) . Representative values for several 
natural surfaces are presented in table 2. These parameters for other 
types of natural surfaces such as brush types, orchard groves, forested 
areas, and urban structure classes are presently unknown. At the present 
time these parameters are believed to be constants for any given surface. 
Throughout this study two natural surfaces of 1 cm lawn grass and 50 cm 
thick grass are compared for the three cases of stability. The very 

u.z 

• * O 

smooth surface of table 2 was not considered since values of >2.5 

V 

00 

are representative of aerodynamic ally rough surfaces. On this basis 
growth of the heated region and decay of surface temperature is presented 
for an arbitrarily selected wind velocity and heat output. 

The superadiabatlc condition extends the heated zone to the great- 
est height and consequently has a more rapid decay of temperature. Thus 
for a given combustion process maximum vertical temperature penetration 
would occur for this condition. A subadiabatic condition experiences 
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a very shallow heated region with a slow temperature decay. This 
condition would experience a more rapid drying out of the surface fuel 
downstream from the combustion zone since the heat generated would be 
concentrated near the surface. 

As the s\irface roughness is increased the importance of stability 
is decreased. The heated zones and temperature patterns tend to approach 
each other. The surface parameters were plotted in figure 15. Extrapo- 
lating to aerodynamic ally rougher sxirface, values of u^u = 0.22 and 
u^z /v = 10^ were selected. The heated zone was compared for a wind 
velocity of 17 ft/sec and is presented in table 5. 


Table 5. -- Growth of heated zone 
uVu =0.22 u„z /v = 10^ u = 17 ft/sec 

o * O 00 00 ' ' 


X 

8 Ft 

Ft 

P = 0.75 

p = 1.00 

P = 1.25 

50 

31.4 

34.8 

24.6 

100 

54.5 

59.0 

60.0 

200 

95.2 

128.0 

142.0 

300 

130.0 

187.0 

260.0 

400 

I64.O 

256.0 

390.0 


Assuming the validity of the extrapolation, atmospheric stability does 
not markedly Influence the growth of the heated region for very rough 
surfaces. 

Influence of stability and siirface roughness on the aerodynamic 
drag is of considerable interest. Measurements of the shear stress at 
the surface have been made by Sheppard (6) and Vehrencamp (8) . Theory 
presented in this study shows that for neutral stability the drag coef- 
ficient is constant. For cases of stability parameter other than unity 
the drag coefficient is a function of the natiiral surface and Reynolds 
Number. The shear stress as a function of wind velocity is shown in 
figxire 14. The case P = 1 is independent of surface type and position. 
For cases of p = O.75 and 1.25 two surface types are con^jared for an 
arbitrarily selected position of 100 ft. For the lower values of veloc- 
ity ratio and macrovlscoslty parameter the shear stress is greatly 
influenced. For the aerodynamically rougher surface the effect of 
stability is decreased. This is consistent with the previous results 
regarding the growth of the heated zone and decay of surface temperature. 

Shear stress results of Vehrencan^) which were obtained on a smooth 
dry-lake bed are presented in figure 14. Theory yields the drag coeffi- 
cient for neutral stability as O.OO36 which is higher by a factor of 3 
when compared with 0.00125 calculated from shear- stress and velocity data 
by Vehrencamp. The shear stress is velocity dependent as specified by 
the stability parameter. This dependency is 

2 

T ~ u 2-p 

O 00 
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VELOCITY RATIO (DIMENSIONLESS) 



Figiire 15. --Relationship between velocity ratio and macroviscosity parameter 


which Indicates the importance of accurate determination of the stability 
of the atmosphere when undertaking field experiments. 

The main step in the theoretical study was the selection of the 
eddy diffusivity for momentum and heat applicable to transport processes 
within the surface air layer. If one assumes the eddy diffusivity for 
momentum and heat as equal and postulates this quantity to be propor- 
tional to the macroviscosity parameter and the normalized distance from 
the siirface to a power, then 



The foregoing expression is a generalized eddy diffusivity which is a 
function of the siirface roughness, distance from the surface, and stabil- 
ity. Experimental verification has been presented by Deacon (2) and is 
here considered the best correlation available at present. Deacon him- 
self, however, expressed some doubt about the applicability of his 
results to correlations of extreme stability or instability. This pro- 
cediire is a step commonly employed by the theorist to produce an accept- 
able mathematical theory. It may be recalled that choosing a constant 
diffusion coefficient (laminar flow case) produces an unacceptable 
theory. 


Certain interesting consequences result from theory. Since q 
is the independent variable given by equation (25), the velocity pro- 
files are dependent on the type of natural siirface, Reynolds Number, 
and stability parameter. The independent variable represents the dimen- 
sionless height and reveals the generalized correlation parameter for 
wind velocity in the atmosphere. As evidenced by figure 4 where 
df/dq = u/u , correlation of the velocity profiles are not universal but 
depend on tSe stability parameter. Likewise, temperature profiles down- 
stream from a heat source in the atmosphere may be correlated by equation 
( 41 ) in terms of the independent variable p . Again figure 5 shows that 
the dimensionless temperatvire profiles are not universal. 

CONCLUSIONS 


1. The kinetic energy rate per unit area of the existing velocity 
field must be much greater than the energy rate per unit area of the com- 
bustion process for transmission of heat by forced convection. P Theory 
yields the condition as 




pu3 

CXD 00 



2. Growth of the heated zone and consequently decay of the temper- 
ature profiles sire markedly affected by atmospheric stability for small 
values of the velocity ratio and macroviscosity parameter. For very rough 
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surfaces, the heated zones and temperature profiles tend to merge for 
all cases of stability. The effect of atmospheric stability on forced 
convection heat transfer over very rough surfaces is unimportant. 


3. The drag coefficient for neutral stability (p = 1) is inde- 
endent of velocity ratio, macroviscosity parameter, and Reynolds Number, 
Theory yields the drag coefficient as 


<^0 = 


= 0.0036 


2 ^ 00^00 


For cases of stability unequal to unity, the drag coefficient is depend- 
ent upon velocity ratio, macroviscosity parameter, and Reynolds Number. 


4. The shear stress at the surface is proportional to the wind 
velocity, l^ccording to 

2 


u 


2-p 


Atmospheric stability does influence the shear stress for small values 
of velocity ratio and macroviscosity parameter but is not very important 
for very rough surfaces. 
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NOMENCLATURE 


The following nomenclature is used in the paper: 


A = constant 


B = constant 

c = specific heat at constant pressure, BTU/lh °F 


Cp = drag coefficient, dimensionless 

e = combustion efficiency, dimensionless 

fCri) = function of independent variable defined by eq. (25), 
dimensionless 

F = function of transformed independent variable, dimensionless 

2 

g = gravitational constant, ft/sec 

g(Ti) = function of independent variable defined by eq. (26), 
dimensionless 

G(ii) = function of independent variable, dimensionless 
H = heating value of fuel, BTU/lb. 

J = mechanical heat equivalent, ft Ib/BTU 
k = Karman constant, O. 4 O 

t o 

k = thermal conductivity of air, BTU/sec sq ft F/ft 

K = constant, equation (32) 

K, = constant, equation (35) 

K(z) = eddy diffuslvity for heat, ft^/sec 



2 

diffuslvity for momentum, ft /sec 


L = line source length, ft 
L* = flame front width, ft 

2 

M = fuel evaporation rate, Ib/ft sec 


N = transformed independent variable, dimensionless 


p = pressure in source zone, psf 
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= atmospheric pressure, psf 
P = dimensionless pressure 
Q = heat output of source, BTU/sec 
R = forest fuel rate, tons/acre min. 

T = temperature in source zone, °R 
= ambient temperature, °R 
= surface temperature, °R 

u = velocity component in horizontal direction, fps 

u = wind velocity, fps 
00 

u^ = friction velocity, fps 
U = dimensionless velocity 

w = velocity component in vertical direction, fps 

W = dimensionless velocity 

X = horizontal distance from source, ft 

x^ = reference length, ft 

X = dimensionless distance 

z = vertical distance above surface, ft 

z^ = roughness length, ft 

Z = dimensionless distance 

Re = Reynolds Number, dimensionless 

Ri = Richardson Number, dimensionless 

P = stability parameter, dimensionless 

5 = height of heated zone, ft 

T) = independent variable, dimensionless 

Tig = independent variable at edge of heated zone, dimensionless 

2 

li = absolute viscosity, lb sec/ft 
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V = ambient kinematic viscosity, ft /sec 
00 

2 , 4 

p = mass density of soiorce zone, lb sec /ft 

2 4 

= ambient mass density, lb sec /ft^ 
a = Prandtl Number, dimensionless 

2 

T = shear stress in source zone, lb/ ft 

2 

T = surface shear stress, Ib/ft 
o 

2 

\|f = stream function, ft /sec 
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